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Ultracold ground state dipolar '*''K*^Rb molecules have recently been produced in a loose harmonic 
trap by employing a magnetic field sweep across a Feshbach resonance followed by stimulated Raman 
adiabatic passage [K.-K. Ni et al, Science 322, 231 (2008)]. The overall experimental efficiency for 
molecule formation was around 20%. We show that the efficiency can be increased to nearly 100% 
if one first loads the atomic gases into an optical lattice of the appropriate depth and tunes the 
interspecies attraction to have exactly one atom of each species at an occupied lattice site. Our 
proposed scheme provides a large enhancement to the dipolar molecule density even at relatively 
high temperatures, and avoids three-body recombination loss by preventing lattice sites from being 
triply occupied. 

PACS numbers: 03.75.Ss 67.85. Pq 67.85.d 71.10.Fd 
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In recent years there has been intense activity in the 
field of ultracold molecules With their long range, 
directional interactions and periodic structure, dipolar 
molecules in optical lattices have attracted particular at- 
tention. Applications of this system include the imple- 
mentation of novel condensed matter phases , quantum 
computation Q and simulation schemes, and assist- 
ing the creation of a dipolar superfluid Q . Motivated by 
these possibilities, a number of recent experiments have 
studied molecules of various isotopic combinations of K 
and Rb [1, 0, i, i], leading to the creation of "OK^^Rb 
molecules in their absolute ground state [l^l. It is de- 
sirable to maximize the efficiency of the association pro- 
cess. When forming heteronuclear molecules from loosely 
trapped gases, phase-space density [llj and the spatial 
overlap of the two components of a Bose-Fermi mix- 
ture have been limiting factors. In an optical lattice, 
molecules may be formed with almost unit efficiency from 
sites containing two atoms, and held for long times with- 
out the possibility of three-body losses [l^. However, 
the limit on the number of molecules that may be made 
with a lattice is the efficiency with which a gas may be 
loaded into it with exactly two atoms per site. 

In this Letter, we show that a two-component gas 
can be loaded into an optical lattice within a harmonic 
trap in such a way that each occupied lattice site con- 
tains one atom of each species. As well as maximizing 
molecule production efficiency, our technique will enable 



the preparation of a lattice with comparatively few de- 
fects, if the experimental entropy per particle can be 
brought to be on the order of the Boltzmann constant, 
kB- This would allow the formation of dense quantum 
degenerate gases of dipolar molecules, and the study of 
interesting new quantum effects where one has both de- 
generacy and long range, directionally dependent inter- 
actions. Such systems could display supersolid phases 
as well as topological insulating behavior. The sys- 
tem we study here is that of fermionic ^"K and bosonic 
^^Rb, which was recently used to produce ultracold dipo- 
lar molecules in their absolute ground state Our 
method, however, is applicable to any two-component 
gas with a comparatively heavy bosonic species and light 
fermionic species. We go beyond previous studies ^ by 
including a realistic trapping potential and the effects of 
finite temperature. 

The experiments of Ref. [l^l took place in a loose har- 
monic trap at temperatures close to the Bose-Einstein 
condensation temperature of *^Rb, Tc- The magnetic 
field sweep to produce Feshbach molecules can have 
near unit efficiency for sufficiently high phase space den- 
sity [l3|, and stimulated Raman adiabatic passage [l3| 
(STIRAP) from the Feshbach molecular state to the ab- 
solute ground state has been measured to have an ef- 
ficiency of 92% [i3, limited by laser phase coherence. 
However, the current efficiency for molecule formation in 
this system (^20%) is limited by the initial phase space 
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TABLE I: Parameters of the Fermi-Bose Falicov-Kimball 
model for a K-Rb mixture and different lattice depths. Note 
that the Rb hopping is more than an order of magnitude less 
than the K hopping and is ignored. The lattice site i lies 
at {aix,aiy), with a the lattice constant, and ix and iy inte- 
gers. The trapping potential y^^'^^ corresponds to a trapping 
frequency of approximately 55 Hz for the 15E^^ lattice, and 
40 Hz for the 20^1?'' lattice. 
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density that can be attained and by the spatial overlap 
of the K and Rb clouds Q . The latter challenge arises 
because K and Rb have different masses but nearly iden- 
tical polarizabilities, so the combined optical and gravi- 
tational potentials yield a lower trap center for Rb. In 
addition, while the Rb cloud shrinks as it is cooled below 
Tc, the K cloud remains large due to the Fermi degen- 
eracy pressure. This further reduces the overlap of the 
clouds, and also limits the ability to sympathetically cool 
the fermions. 

One possible approach to achieving a high joint den- 
sity of K and Rb would be to simply load the gases into 
a deep optical lattice, so that the Rb atoms form a Mott 
insulator and the K atoms form a spin-polarized band 
insulator. This approach would ensure uniform density 
of both species at the center of the trap. However, a 
substantial fraction of both species would remain in the 
outer region of the trap, where K and Rb have "metal- 
lic" and superfluid character, respectively, which produce 
a nonuniform distribution. Furthermore, the sizes of the 
K and Rb insulating regions will likely not match. Con- 
sequently, this approach would not be able to pair all 
K and Rb atoms, and would leave a severe limitation on 
the possible efficiency of molecule formation. Instead, we 
propose to load the atoms into an optical lattice with the 
Rb atoms deep in the Mott phase but the K atoms able 
to hop between sites, and tune the interspecies attraction 
into a range where the energetically stable phases have 
precisely one K and one Rb atom per lattice site. In this 
case, at low temperature, the system will minimize its en- 
ergy by creating these pre-formed molecules, rather than 
having only K or Rb atoms on a lattice site. The system 
will then be perfectly primed for the formation of ground 
state molecules with a Feshbach sweep and STIRAP. 

We consider optical lattices that are deep enough for 
the tunnelling of Rb atoms to be negligible. We may 
then neglect the kinetic energy and describe the spin- 
polarized K-Rb mixture with the Fermi-Bose Falicov- 



Kimball model [H, H [13 , 

ij i 

+ u^~^'Y.4cMb. (1) 
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where cj creates a fermionic ^"K atom at lattice site i and 
creates a bosonic ^^Rb atom at site i. The parameters 
include the K hopping i^, trap potential (V^^) and chem- 
ical potential (/i^), as well as the corresponding terms for 
the Rb atoms. In addition, we have the attractive inter- 
species interaction strength, U^~^^, and the repulsive 
Rb-Rb interaction strength, jjRb-Rb ^ The Rb scattering 
length is lOOao, where ao = 0.0529nm The K-Rb 

scattering length varies with magnetic field strength B 
as -191ao(l - ^B/[B - B^]), where Bq = 546.9 G and 
AB = — 3.1 G are the location and width of the resonance 
employed to make the Feshbach molecules, respectively, 
and IG = 10-4 T 0. The desired U^'^'' is set by 
choosing the appropriate value of B close to Bq. 

We assume a 1030 nm laser for the optical lattice, 
which is chosen because it provides approximately the 
same trapping potential for K and Rb [6|. We work on 
a two-dimensional lattice with 51 x 51 lattice sites. The 
particle number is fixed to be 625 K atoms and 625 Rb 
atoms in the two-dimensional plane, which is close to half 
filling of the simulated lattice. Wc examine two lattice 
depths which are both deep in the Mott insulating phase 
for Rb. The first is 15£'|'' and the second is 20£:|J^ 
where ~ /i x 2.164 kHz is the recoil energy of the Rb 
atoms. Assuming this lattice is created by introducing a 
AQE^ potential in the third dimension, the parameters 
of the Hamiltonian are summarized in Table |l] for the two 
lattice depths. 

The phenomenon of pre-formed molecules in the 
Falicov-Kimball model has been widely studied for the 
case of hard-core bosons (or equivaleiitly spin polarized 
fermions) on a homogeneous lattice [l9|. In addition, 
the case of fermions in a harmonic trap was examined 
in Ref. Here we extend this work to examine soft- 

core bosons in a harmonic trap. We focus on two spe- 
cific quantities. The efficiency of producing pre-formed 
molecules measures the joint probability to find exactly 
one K and exactly one Rb atom on a lattice site. The 
entropy is employed as an effective temperature scale, 
under the assumption that the optical lattice is turned 
on adiabatically and that the lattice does not cause sig- 
nificant heating during the experiment. We have two 
complementary calculational tools that we employ. Inho- 
mogeneous dynamical mean-field theory (IDMFT) ,21;] is 
approximate for a two-dimensional system, but allows us 
to determine both the production efficiency of pre-formed 
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FIG. 1: (Color online.) Efficiency for pre-forming molecules 
on the 15i?|{'' lattice for two different [J^-^^ interactions. 
Filled symbols are IDMFT and open symbols are QMC. 
^K~Rb ^ _2^K _iot^' correspond to B = 542.9 G and 
565.9 G, respectively. 



molecules and the entropy per particle of the mixture. 
In contrast, direct quantum Monte Carlo (QMC) simula- 
tions (2^ are numerically exact, subject to statistics and 
the possibilities of becoming stuck in local minima, but 
can only accurately determine the efficiency. We em- 
ployed a zero-temperature local-density-approximation 
analysis to determine the set of parameters that we stud- 
ied at finite temperatures. 

In Fig. [1] we plot the efficiency versus temperature for 
the two different calculational methods on the 15E^'' lat- 
tice. The horizontal axis is on a logarithmic scale to 
emphasize the differences between the two curves, which 
disagree only in the narrow region where the efficiency 
changes rapidly. We note that the disagreement between 
the two approaches is smaller for deeper lattices. Con- 
sequently, the results of Fig. [1] verify that the IDMFT 
gives accurate results in the parameter range considered 
here. Another consistency check we used was to exam- 
ine the sum rules for the first three moments of the local 
density of states and for the zeroth and first moment of 
the imaginary part of the self-energy |23| . The sum rules 
are satisfied to high accuracy for all cases except very 
low temperature, where these functions have sharp peaks 
that are not easily represented on our coarse frequency 
grid. This further confirms the quality of the IDMFT 
approximation. 

In Fig. [5] we report our main results for efficiency ver- 
sus entropy per particle at the two different optical lat- 
tice depths that we consider. The entropies shown corre- 
spond to temperatures ranging from 0.05t^ to lOt^. In 
each case we consider three different values of the inter- 
species interaction strength. All of the values used are 
weak enough that the number of sites with two Rb atoms. 
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FIG. 2: (Color online.) Efficiency for pre-forming molecules 
on the (a) 15E^'' and (b) 20E^'' lattices versus the entropy 
per particle, as calculated with IDMFT. The different sym- 
bols correspond to different K-Rb interaction strengths. The 
efficiency is nearly 100% at a sufficiently low entropy per par- 
ticle, corresponding to low temperature, and for large enough 
K-Rb attraction. 



and accompanying three-body recombination loss, are 
acceptably low. We note that for a fixed value of the 
entropy per particle, the deeper lattice will provide a 
higher efficiency. For a given lattice depth, once the 
attraction is strong enough to stabilize the prc-formed 
molecule, the efficiency depends only weakly on the in- 
teraction strength. By tuning the scattering length into 
this regime, one can increase the efficiency to 90% in the 
shallower lattice, and to more than 95% in the deeper 
lattice, even for the relatively high temperature corre- 
sponding to an entropy of fcs per particle. 

In Fig. [3] we plot the densities of the Rb and K atoms 
for the deeper lattice with jj^^Rb = —IQt^. The size 
of the symbols indicate the density on each lattice site. 
Note how the Rb and K atoms generally lie on top of 
each other, but the density approaches one only at the 
lower temperatures. At higher temperatures, it is more 
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FIG. 3: (Color online) Density distributions for Rb (left) 
and K (right) atoms, at temperatures indicated on the plots 
along with the resulting efflciency of producing preformed 
molecules. Here, the lattice has a depth of 2QE^ , and the 



interspecies attraction is U 



= — 16t . The symbol size 



is proportional to the density on each lattice site. Sites inside 
the solid circle, at the center of the top left panel, show the 
lattice sites that have a double occupancy of Rb atoms more 
than 2% of the time. 



likely that a site could have two Rb atoms. We indicate 
the temperature at which this begins to become notice- 
able with the solid circle in Fig. [3l within which sites are 
doubly occupied more than 2% of the time. At the low- 
est temperatures, we can see how the reduced density of 
the cloud in the outer regions corresponds to nontrivial 
ordered patterns of the pre-formed molecules, i.e. that 
the few occupied sites still have one K atom and one 
Rb atom. One can also see that as the temperature is 
made too large, there are boundary effects due to the fi- 
nite size of our system, which restricts the clouds from 
becoming as large as they should be. This effect overes- 
timates the efficiency and corresponds to the rightmost 
data in Figs. [1] and [H 

In conclusion, the production efficiency of ground 
state dipolar molecules can be greatly increased by first 



loading the atomic gases into an optical lattice with 
an interspecies scattering length such that the system 
tends to form sites containing one atom of each species. 
Our scheme works well even for gases prepared at rela- 
tively high temperatures, but is predicted to improve as 
this temperature decreases. Finding the optimal lattice 
depth, however, requires a number of different choices to 
be made. Although our calculations predict higher ef- 
ficiency for deeper lattices, experimental issues such as 
the longer times required for adiabatic loading require 
consideration. As the range of scattering lengths within 
which pre-formed molecules are stable varies with the lat- 
tice depth, optimizing these parameters requires further 
input from experiment. 
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